Effects of nanoscale spatial inhomogeneity in strongly correlated systems 
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We calculate ground-state energies and density distributions of Hubbard superlattices character- 
ized by periodic modulations of the on-site interaction and the on-site potential. Both density-matrix 
renormalization group and density-functional methods are employed and compared. We find that 
small variations in the on-site potential Vi can simulate, cancel, or even overcompensate effects due 
to much larger variations in the on-site interaction Ui. Our findings highlight the importance of 
nanoscale spatial inhomogeneity in strongly correlated systems, and call for reexamination of model 
calculations assuming spatial homogeneity. 
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A large part of the complexity of strongly correlated 
systems arises from the multiple phases that coexist or 
compete in their phase diagrams. Metallic and insulat- 
ing phases are separated by metal-insulator transitions, 
and subject to the formation of various types of long- 
range order, such as antiferromagnetism, superconduc- 
tivity, and charge or spin-density waves. The relative 
stability of such phases is determined by differences in 
appropriate thermodynamic potentials, or, at zero tem- 
perature, in their ground-state energies. Identification of 
the appropriate order parameters and calculation of the 
ground-state energies of the various phases is a compli- 
cated problem, and the nature of the phase diagram of 
many strongly correlated systems is still subject to con- 
siderable controversy. It is widely believed, however, that 
a minimal model containing the essence of strong cor- 
relations, and displaying many of the above-mentioned 
phases, is the homogeneous Hubbard model, which in 
one dimension and standard notation reads 

H hom = -ty^(ct ff ^ + i, -+-g.c.) + t/y^ct t Ct-|-ct | Cij p . (1) 

i,a i 

Much theoretical effort is thus going into the analysis of 
the homogeneous Hubbard model and the clarification of 
the nature of its ground state. 

In a parallel development, nanoscale spatial inhomo- 
geneity has been observed experimentally to be a ubiqui- 
tious feature of strongly correlated systems jiS^ii^LSi 
but although its importance is widely recognized, the 
consequences of such inhomogeneity are still insufficiently 
understood. The present paper investigates the effects 
of, and the competition between, two different manifes- 
tations of nanoscale inhomogeneity in strongly correlated 
systems: local variations in the on-site potential and in 
the on-site interaction. We base our analysis on the in- 



homogeneous Hubbard model 

Hinhom = — t'5^X C ia Ci + 1 <<r + H.C.) 

+ Yl U A^\ C \i^i + V i C iv C i<r, (2) 
i ia 

which differs from the homogeneous model Q by allow- 
ing for spatial variations in the on-site interaction Ui 
and the presence of the on-site potential Vi. Variations 
in Ui and Vi may arise, e.g., due to inequivalent sites 
in the natural unit cell, modulation of system parame- 
ters in artificial heterostructures, or self-consistent mod- 
ulations in local system properties due to formation of 
charge-ordered states. In this paper we are specifically 
concerned with one-dimensional superlattice structures 
in which both Ui and Vi vary periodically on a length 
scale comparable to, or somewhat larger than, the lat- 
tice constant. Such superlattices have recently attracted 
much attention due to their complex ground-state and 
transport properties. 9 ! 10 ! 11 ' 12 : 13 ! 14 : 15 ! 16 : 17 ! 18 Our results, 
reported below, have a direct bearing on the investigation 
of such superlattices. However, for our present purposes 
the most important aspect of superlattice structures is 
that they constitute a representative system in which the 
consequences of nanoscale spatial variations of system pa- 
rameters in the presence of strong Coulomb correlations 
can be explored systematically. Accordingly, we expect 
our main conclusions to hold also in many other spatially 
inhomogeneous correlated systems. 

Figure^shows the density profile of a typical superlat- 
tice structure in which the on-site interaction U is mod- 
ulated in a repeated pattern of repulsive (U — 3) and 
nonintcracting (Ui — 0) 'layers' with Ljj and Lq sites, 
respectively, and the on-site potential Vi is taken to be 
constant at all sites. The two curves shown were obtained 
with different many-body techniques. The full curve 
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TABLE I: Ground-state energy of two open superlattices 
with modulated on-site interaction Ui and spatially constant 
on-site potential Vi, obtained with DMRG and with DFT/BA- 
LDA. Upper part: large lattice with L — 300 sites. Ly = 10 
interacting sites (Ui — 3) alternate with Lo = 10 noninter- 
acting sites (Ui = 0). Lower part: strongly modulated lattice 
with L — 100 sites. Lu — 1 interacting site (Ui — 6) al- 
ternates with Lo = 1 noninteracting site. N is the number 
of fermions, and the column labeled A% contains the abso- 
lute percentual deviation of the DMRG from the DFT/BA- 
LDA values. The agreement between DMRG and BA-LDA is 
slightly better for the more slowly modulated lattice. 25 
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was obtained using the density-matrix renormalization 
group (DMRG) while the dotted curve was obtained 
from density-functional theory (DFT) within the Bcthc- 
Ansatz local-density approximation (BA-LDA) ,21*22^ In 
view of the complexity of the problem and the surprising 
nature of some of our conclusions, we found it advisable 
to bring two independently developed and implemented 
many-body methods to bear on the problem. 

DMRG is a well-established numerical technique, 
whose precision can be improved systematically, at the 
expense of increased computational effort In our 
DMRG calculations, truncation errors were kept of the 
order of 10~ 6 or smaller, and increasing the precision 
beyond this did not affect any of our conclusions. BA- 
LDA is a more recent developmenfeSLSSiS^ (although the 
original LDA concept is, of course, widely used in ab ini- 
tio calculations). In LDA calculations the final precision 
is ultimately limited by the locality assumption inher- 
ent in the LDA, and improvements must come from the 
development of better functionals. This intrinsic limita- 
tion of LDA is offset by its applicability to very large 
and inhomogeneous systems, at much reduced compu- 
tational effort: Calculations for the type of superlattice 
structures investigated here typically take only seconds to 
minutes with BA-LDA, regardless of the type of bound- 
ary condition usedi^i Final BA-LDA results for densities 
and energies typically agree with DMRG ones to within 
< 3%, the agreement being slightly better for energies 
than for densities^ Here we consider both methods as 
complementary. All essential conclusions reported below 
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FIG. 1: Density profile of a one-dimensional superlattice with 
L — 60 sites, N = 30 fermions, open boundary conditions, 
and a superlattice structure consisting of a periodic sequence 
of Lu = 3 interacting (Ui = 3) and Lo = 2 noninteracting 
(Ui = 0) sites. Full curve: DMRG calculation. Dotted curve: 
DFT/BA-LDA calculation. 

were obtained on the basis of independently implemented 
and performed BA-LDA and DMRG calculations. As an 
illustration, Table [3 compares ground-state energies ob- 
tained with both methods for one much larger and one 
much more rapidly modulated superlattice than the one 
shown in Fig. ^ 

Inspection of Fig. [3 shows that an attractive potential 
on the repulsively interacting sites can completely reverse 
the effect of the Coulomb repulsion Ui and draw a sub- 
stantial number of electrons to the interacting sites (cir- 
cles in Fig. |2J). While this might have been anticipated 
qualitatively as a result of the competition between an 
attraction and a repulsion, it comes as a surprise that 
the effect of the (often neglected) variations in the on- 
site potential is much stronger than the one of variations 
in the on-site interaction: already a very weak attrac- 
tive potential suffices to smooth out the density distri- 
bution, resulting in an essentially homogeneous charge 
profile (triangles in Fig. 0). Although we have taken su- 
perlattices as our example, the effect is clearly not de- 
pendent on periodicity of the modulations in Ui and Vi, 
and is expected to show up rather generally. 

Figure ^ an d Table |U represent superlattices in 
which only the on-site interaction is spatially mod- 
ulated, which is the case mostly studied in the 
literaturei^ 9 ^ ! 1 ^ 12 ! 13 ^ 4 ! 15 ^ 6 ! 17 Many important features 
of superlattice structures are already apparent in this 
type of model. However, in a real system it is impossible 
to modulate the on-site interaction without simultane- 
ously modulating the on-site potential as well, i.e., with- 
out creating inequivalent sites. Such a double modula- 
tion is found, e.g., in artificially grown layered structures, 
in impurity systems, and in periodic arrays of Fermi- 
liquid leads (corresponding to approximately noninter- 
acting sites) and quantum wires/dots (corresponding to 
interacting sites). Fig. [21 illustrates the consequences a 
modulation of the on-site potential has on the density 
profile of a system in which both interaction and poten- 
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FIG. 2: Density profiles, obtained with BA-LDA^ of an 
L = 100 site system of N — 50 fermions with periodically 
modulated on-site interaction of amplitude U = 3 (Lu = 
I/O = 10), and periodic boundary conditions. Squares: no 
modulation in on-site potential (vi = 0). Circles: on-site 
potential modulated such that Vi = —2 on the interacting 
sites and Vi = on the noninteracting sites. The density 
profile is inverted, indicating overcompensation of Ui by Vi. 
Triangles: on-site potential modulated such that Vi = —0.555 
on the interacting sites and Vi — on the noninteracting sites. 
The superlattice structure is erased from the density profile. 
The lines are guides for the eye. 



tial vary. 

In Fig. |3 we compare the Friedel oscillations arising 
from the system boundaries in a homogeneous system 
with the ones arising in a superlattice of same size and 
with the same number of fermions, but subject to pe- 
riodic modulations of Ui and Vi, chosen such that both 
density profiles become similar. We have deliberately not 
chosen modulation parameters that optimize the agree- 
ment between both curves, because had we done so they 
would be visually indistinguishable on this scale. 

Figure |2 shows that a small value of Uj can have 
stronger effects than a larger value of E7j, while Fig. 
shows that essentially the same modulation pattern in the 
density profile can be obtained from either v or U. These 
are unexpected findings. Normally it is assumed that 
in systems modeled by the Hubbard model the particle- 
particle interaction U is much more important than the 
on-site potential, which is mostly taken to be spatially 
constant, or, if it varies, to produce only minor modifi- 
cations in systems whose physics is governed by U. 

To investigate in more detail this competition between 
on-site interaction and on-site potential we need to es- 
tablish criteria for comparing the consequences of Vi and 
of Ui. Motivated by the experimental observation of 
nanoscale density variations and by the theoretical im- 
portance of ground-state energies for analyses of phase 
diagrams, we adopt two distinct criteria. Criterium (i) 
consists in searching for that modulation of the on-site 
potential vi in a doubly modulated lattice that cancels 




FIG. 3: Full curve: density profile, obtained with BA-LDA,2 6 
of a L = 160 site homogeneous system with open boundary 
conditions, Ui = 2, N — 80, Vi — 0. Friedel oscillations aris- 
ing from the system boundaries are clearly visible. Dotted 
curve: density profile of same system but subject to modula- 
tions periodically alternating Lu = 6 sites with Ui = 2 and 
Vi — —0.2, with La = 10 sites with Ui — Vi — 0. The super- 
lattice structure due to the presence of both modulations is 
completely erased from the density profile, while the Friedel 
oscillations arising from the boundary remain prominent. 



the effect of the modulation of Ui on the density distribu- 
tion, i.e., smoothes out the oscillations, making the net 
density homogeneous. A particular example of this can- 
cellation is given by the triangles in Fig. [5] Criterium (ii) 
consists in searching for that modulation of the on-site 
potential Vi in the doubly modulated lattice that yields 
the same ground-state energy Eq as in a homogeneous 
lattice with Vi = and Ui = U at all sitesiSi Our results, 
displayed in Fig. 01 24 show that, regardless of whether one 
adopts the density or the energy criterium, the modula- 
tion of the on-site potential required to cancel the effect of 
the modulation of the on-site interaction is up to an order 
of magnitude smaller than U. For open boundary con- 
ditions we have obtained the same conclusion also from 
DMRG calculations. Changes in the modulation pattern 
do not change the order of magnitude of the ratio of \v\ 
to U appreciably. 

A semi-quantitative explanation for this relation of \v\ 
to U can be given within DFT, by considering the effec- 
tive potential entering the Kohn-Sham equations for the 
Hubbard model, u e //,i — Vext.i +vn.i + v c i . For unpolar- 
ized systems (n-\^ — n^i — rij/2) the Hartree potential 
VH,i can be written v H ,i = Uiuijl. Within BA-LDA DFT 
the density and total energy are thus calculated from 
an effective Hamiltonian containing the modulated inter- 
action and external potential only via the combination 
v e xt,i + Uirii/2 + v Ct i(ni, Ui). Since the correlation poten- 
tial v c< i is typically about an order of magnitude smaller 
than Vn,i, the modulated interaction Ui enters the effec- 
tive Hamiltonian approximately on the same footing as 
the modulated potential v ex t,i, but renormalized by the 
factor ni/2. In the above calculations the average density 
N/L = 0.5, and the local density n, is not very differ- 
ent. The upshot is that self-consistent screening of the 
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FIG. 4: Open squares: Amplitude of the modulation of the 
attractive on-site potential that cancels as much as possible 
the effect of the modulation of the repulsive on-site interac- 
tion, leading to a homogeneous density profile [criterium (i)]. 
Full circles: Amplitude of the modulation of the attractive on- 
site potential that reproduces in the doubly modulated system 
the energy of the homogeneous system [criterium (ii)]i— Sys- 
tem parameters: L — 160 sites, TV = 80 fermions, Ljj = 6 in- 
teracting sites, alternating with Lo = 10 noninteracting sites. 



particle-particle interaction effectively reduces the mod- 
ulation in the interaction by a factor ~ 4, compared to 
modulations in the potential, in good agreement with the 
numerical results in Fig. 

Of course, in a real system one cannot adjust t>, at will, 
and the precise fine-tuning required to obtain smooth 
density profiles, or energies identical to the ones found 
in homogeneous systems, is not expected to occur fre- 
quently in nature. The main implication of these criteria 
is rather that they establish a scale for comparison of v 



and U, indicating that even weak spatial variations of v 
can be more important than much stronger ones in U . 
This observation flags a warning signal to the use of ho- 
mogeneous Hubbard models (or ones in which only U is 
modulated) in the analysis of situations characterized by 
nanoscale spatial inhomogeneity, such as the pseudo-gap 
phase of cupratesiSiiMiiiiliiL 8 . or superlattices and similar 
heterostructures, 9 ! 10 ! 11 ! 12 . 13 ! 14 ! 15 ! 16 ! 17 ! 18 

We conclude that even in the presence of strong corre- 
lations, spatial variations of the on-site potential Vi are 
not a minor complication in a system dominated by the 
on-site interaction Ui, but a major effect, which crucially 
contributes to observables, and can mask or overcom- 
pensate the effect of the interaction on the density pro- 
file, ground-state energy, and other quantities. For the 
density profile, this means that attempts to model the 
microscopically inhomogencous charge distribution, seen 
experimentally, 1 ' 2 ' 3 !^*^* 7 * 8 , by Hubbard models that are 
homogeneous or that modulate only the interaction Ui, 
cannot lead to conclusive results. The influence of mod- 
ulations in Vi on the ground-state energy, on the other 
hand, implies that an analysis of the relative energetic 
stability of the various phases appearing in strongly- 
correlated systems is incomplete, and potentially mis- 
leading, if the effects of spatial inhomogeneity in these 
phases are not taken into account. All this calls in ques- 
tion the common practice to employ the homogeneous 
Hubbard model to model spatially inhomogeneous many- 
body systems, and demands a reconsideration of the role 
of nanoscale spatial inhomogeneity in strongly correlated 

S y S t eIIls 1-2. 3.4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18 

Acknowledgments 

This work was supported by FAPESP and CNPq. 



* Electronic address: capelle@if.se. usp. 

1 'Intrinsic Multiscale Structure and Dynamics in Complex 
Electronic Oxides', A. R. Bishop, S. R. Shenoy and S. Srid- 
har, eds. (World Scientific, New Jersey, 2003). 

2 K. McElroy, D.-H. Lee, J. E. Hoffman, K. M. Lang, 
E. W. Hudson, H. Eisaki, S. Uchida, J. Lee, and 
J. C. Davis, cond-mat/0404005 (2004). 

3 D. J. Derro, E. W. Hudson, K. M. Lang, S. H. Pan, 
J. C. Davis, J. T. Markert, and A. L. de Lozanne, Phys. 
Rev. Lett. 88, 097002 (2002). 

4 A. Yazdani, C. M. Howald, C. P. Lutz, A. Kapitulnik, and 
D. M. Eigler, Phys. Rev. Lett. 83, 176 (1999). 

5 E. W. Hudson, S. H. Pan, A. K. Gupta, K.-W. Ng, and 
J. C. Davis, Science 285, 88 (1999). 

6 K. M. Lang, V. Madhavan, J. E. Hoffman, E. W. Hudson, 
H. Eisaki, S. Uchida, and J. C. Davis, Nature 415, 412 
(2002). 

7 S. H. Pan, J. P. O'Neal, R. L. Badzey, C. Chamon, H. Ding, 
J. R. Engelbrecht, Z. Wang, H. E. Eisaki, S. Uchida, 
A. K. Gupta, K.-W. Ng, E. W. Hudson. K. M. Lang, and 
J. C. Davies, Nature 413, 282 (2001). 

8 T. Cren, D. Roditchev, W. Sacks, and J.Klein, Euro. Phys. 
Lett. 54, 84 (2001). 



9 A. L. Malvezzi, T. Paiva, and R. R. dos Santos, Phys. 
Rev. B 66, 064430 (2002). [Virtual Journal of Nanoscale 
Science & Technology 6 (2002) (www.vinano.org)]. 

10 J. Silva- Valencia, E. Miranda, and R. R. dos Santos, Phys. 
Rev. B 65, 115115 (2002). 

11 J. Silva- Valencia, E. Miranda, and R. R. dos Santos, J. 
Phys. Condens. Mat. 13, L619 (2001). 

12 T. Paiva and R. R. dos Santos, Phys. Rev. B 65, 153101 
(2002), ibid 62, 7007 (2000), ibid 58, 9607 (1998). 

13 Phys. Rev. Lett. 76, 1126 (1996). 

14 T. Paiva, M. El. Massalami, and R. R. dos Santos, J. Phys. 
Condens. Mat. 15, 7917 (2003). 

15 D. Gora, K. J. Rosciszewski and A. M. Oles, J. Phys. Con- 
dens. Matter 10, 4755 (1998). 

16 M. Noh, D. C. Johnson and G. S. Elliott, Chem. Mater. 
12, 2894 (2000). 

17 L. Bauernfeind, W. Widder, H.F. Braun, Physica C 254, 
151 (1995). 

18 S. G. Ovchinnikov, Phys. Usp. 46, 21 (2003). 

19 S. R. White, Phys. Rev. Lett. 69, 2863 (1992), ibid 77, 
3633 (1993). 

20 A. L. Malvezzi, Braz. J. Phys. 33, 55 (2003). 

21 N. A. Lima, M. F. Silva, L. N. Oliveira, and K. Capelle, 



5 



Phys. Rev. Lett. 90, 146402 (2003). 

N. A. Lima, L. N. Oliveira, and K. Capelle, Europhys. 
Lett. 60, 601 (2002). 

K. Capelle, N. A. Lima, M. F. Silva, and L. N. Oliveira, 
in ' The fundamentals of density matrix and density func- 
tional theory in atoms, molecules, and solids', N. Gi- 
dopoulos and S. Wilson eds. (Kluwer, Dordrecht, 2003) 
cond-mat/0209245 . 

The computational efficiency of BA-LDA is crucial for the 
production of Fig. 2] which requires for each value of Ui 
many calculations with different values of Vi, until the Vi 
that smoothes the density or equalizes the energy is found. 
Improved LDA functionals, currently under construction, 
have in preliminary calculations shown even smaller devi- 
ations from DMRG data. 

For visual clarity, Figs.|5|to0]contain only curves obtained 



with BA-LDA, but as Fig. [Hand Table Q] show, the differ- 
ence between DMRG and DFT/BA-LDA results is much 
smaller than the difference between the effect of Ui and 
that of Vi, which is robust and not dependent on compu- 
tational methodology. 

Note that one cannot use the Hohenberg-Kohn theorem to 
argue that the potential that produces the same density 
distribution in both systems must be the same that yields 
the same ground-state energy. The ground-state energy is a 
unique functional of the density only for fixed interaction, 
whereas in the present case the two superlattices that are 
being compared have different interactions. Hence the two 
criteria are different, and the two sets of potentials in Fig. [I] 
are not identical. 
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